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Abstract 

Codimension two branes play an interesting role in attacking the cosmological constant problem. 

f"^ ' Recently, in order to handle some problems in codimension two branes in Einstein gravity, Bostock 

O ' 

ps) ' et al. have proposed using six-dimensional Einstein-Gauss-Bonnet (EGB) gravity instead of six- 

S , dimensional Einstein gravity. In this paper, we present the solutions of codimension two branes in 

six-dimensional EGB gravity. We show that Einstein's equations take a " factorizable" form for a 
factorized metric tensor ansatz even in the presence of the higher-derivative Gauss-Bonnet term. 

^ ■ Especially, a new feature of the solution is that the deficit angle depends on the brane geometry. 

o. 

C^^ , We discuss the implication of the solution to the cosmological constant problem. We also comment 

^O I on a possible problem of inflation model building on codimension two branes. 
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I. INTRODUCTION 

The idea of braneworlds and large extra dimensions [l| implies that the cosmological 
constant problem (see Refs. [^, y] for reviews) may be a clue of our unawareness of the true 
nature of spacetime: vacuum energy may be large, but it simply does not gravitate in the 4- 
dimensional braneworld where we are living in. The key point is that, cosmological constant 
is a reflection of 4-dimensional spacetime geometry and thus is what is observed directly in 
cosmological observations. The puzzle arises only after we use General Relativity to find 
that cosmological constant describes vacuum energy of standard model particles. So if we 
modify gravity theory by introducing higher dimensional spacetime and objects like branes, 
it is possible that 4-dimensional cosmological constant is not linked to 4-dimensional vacuum 
ene.,, bu. scetKing else such a. Kigh. di.e.ica. vacuu. eue.g. See. e.g.. Refs.Q Q 
for some earlier endeavors in this direction. 

Recently, Carroll and Guica presented an interesting exact solution of this type p. They 
considered a factorizable braneworld spacetime with two extra dimensions and explicit brane 
sources. The compactification manifold has the topology of a two-sphere, and is stabilized 
by both a bulk cosmological constant and a magnetic flux. From their solution, they found 
that the fiat nature of the 4-dimensional geometry is independent of the brane tension. This 
feature moves the cosmological constant problem completely into the extra dimensions. Of 
course, this is not a complete solution to the cosmological constant problem, since it still 
needs fine-tuning in the bulk. But it transforms the nature of the problem in a suggestive 
way. 

The interesting feature of Carroll and Guica's solution is not an accident (see Ref. ItII for 

n 

some other models of codimension 2 branes that share the similar feature, see also Ref. [8|] 
for earlier ideas along this line). It can be shown that the independence of 4-dimensional 
geometry on the brane tension is a general feature of codimension 2 branes in factoriz- 
able spacetime in Einstein gravity. The following discussions will also be helpful for us to 
understand the properties of codimension 2 branes in Einstein-Gauss-Bonnet gravity (see 
SecUH). 

Let's consider a factorizable metric ansatz, 

ds^ = GAEdX^dX"" = g^,{x)dx^dx^ + ^ab{y)dy''dy'' (1) 

where A,B = 0, ...,5, fJ.,i^ = 0, ...,3 and a,b = 4,5. The •jab is the metric of a Einstein 



manifold with curvature k = —1,0,1. Note that due to the presence of branes, there will be 
deficit angles in the extra dimensions at the positions of the branes (see Sec JIII|) . but this 
will not influence the local geometry of the extra dimensions at other points. 

We will consider the simplest model of branes which is also the case considered in most 
of the literature on braneworld cosmology: the branes are described by Nambu-Goto action 
(see Ref.J9| for an elegant review), 



SnG — I d X^/\G\Cbrane , (2) 

where 



Chrane = -J^ J ^^xJ -^a,6^'HX ~ X,{x)) , (3) 



in which i labels the branes, cij and Xj are the tension and position of the ith brane, 
respectively. The energy-momentum tensor of branes follow by varying G"^^ in (j21) 9| 

^^B = -E^f'"°V^'Hy-y.), (4) 

iV^\o oj 

With the help of the fact that the Einstein tensor. Gab = Rab — \habR, vanishes identically 
for any 2-dimensional metric hab and denoting the bulk energy- momentum tensor by T^^, 
contracting the transverse component of the Einstein equation gives 

m = -^Ti , (5) 

where T^ = T^7"^; while contracting the longitudinal component gives 

fiH + iflbl = -5^|rf + T']. (6) 

where T^ = Tj^^g^^ and T^ = T'^^g^^'^ . Now, Eq.(jS)) tells us that the scalar curvature 
of the ^-dimensional spacetime is totally determined by the transverse component of the 
total bulk energy-momentum tensor. Thus if we assume the 4-dimensional geometry to be 
maximal symmetric, then it is determined totally by the transverse part of the bulk energy 
momentum tensor. Specifically, the 4-dimensional geometry does not depend on the brane 
tension. Then, after we find R[g] from Eq.(jSI), substituting it into Eq.®, we can find 
the bulk curvature -R[7]. In sum, for codimension 2 branes in factorizable spacetime, the 
brane geometry is determined by the transverse component of the Einstein equations and 



the bulk geometry is determined by the longitudinal component of the Einstein equations. 
Roughly speaking, we can say that Einstein equations in factorizable spacetime are also 
"factorizable" . This is the secret of codimension 2 branes in Einstein gravity. 

While the above discussion is exiting, unfortunately, when considering realistic cosmolog- 
ical evolution of this model, we will encounter some fundamental difficulties. One of them 
is that if we assume the brane energy-momentum tensor to be of the form as the perfect 
fluid, i.e. Tj^ = {PjP, p,p}5^, then p and p must satisfy p + p = 0, i.e. it behaves like 
the brane tension |lO|. This forbids us adding dust and radiation on the brane, thus it is 
cosmologically unrealistic. To remedy this and other difficulties of codimension 2 branes 
in Einstein gravity, recently, Bostock et al. suggested that we may add the Gauss-Bonnet 



term to the 6-dimensional gravitational action [11[ (however, see also Ref. 1^ and reference 
therein for some other suggestions to handle this problem). It is also worth commenting 
that the idea that 6-dimensional Einstein- Gauss-Bonnet (EGB) gravity could be relevant in 
relation to the cosmological constant problem was originally presented in Ref.jlJ] (see also 
Ref. [lJ| for some subsequent related works). The Gauss-Bonnet term is quadratic in the cur- 
vature tensors and is a topological invariant in 4-dimensional manifold (see, e.g., Ref.(l5l|): 
but in higher dimensions, it has the well-know property that the equation of motion derived 
from it remains second order differential equations of the metric. Furthermore, considering 
higher derivative terms is also necessary to develop the braneworld scenario in a more string 
theoretic setting (see, e.g., Ref. |lfi|). Specifically, the Gauss-Bonnet combination arises as 
the leading order for quantum corrections in the heterotic string effective action and is the 
only quadratic combination of curvature tensors that is ghost-free jnf . 

Thus, the investigation of codimension 2 branes in EGB gravity is well motivated (see 
Ref. [18[ for some other recent discussion of codimension 2 branes in EGB gravity). Of 
course, one of the best way to understand the property of a gravity theory is studying its 
exact solutions. Especially in the present case, the EGB gravity is intended to remedy the 
model in Einstein gravity. Thus one natural step is to derive and compare the corresponding 
solutions in EGB gravity under the same assumption of spacetime geometry and matter 
content with Einstein gravity case. In particular, it is important to check that the important 
property in Einstein gravity, i.e. the independence of the 4-dimensional geometry on the 
brane tension, is retained in EGB gravity. If this were not the case, considering EGB gravity 
would not be so well motivated. We will see in Secln]that the discussion above for Einstein 



gravity also applies to EGB gravity, thus EGB gravity retains the main features of Einstein 
gravity. In Sec lIIIl we will also see that some new features will arise in EGB gravity. The 
last section, Sec ilVl is devoted to conclusions and we comment on inflation model building 
in the codimension 2 brane scenario. 

II. EINSTEIN-GAUSS-BONNET EQUATION IN FACTORIZABLE SPACETIME 

Let's consider adding the Gauss-Bonnet term to modify the 6-dimensional gravity jll|, 
which is described by the action 

Se = J d'X^\^[R + aRU , (7) 

where a is the Gauss-Bonnet coupling constant with dimension [a] = (mass)^'^. Following 
the original derivation [1^, one generally assumes a > 0, but in the literature the a < 
case is also often discussed. We will see in Sec JIIII that, from the exact solution we found, 
the requirement of the geometry to be nonsingular will rule out a negative Gauss-Bonnet 
coupling constant. The Gauss-Bonnet term Rgb is given by 

p2 _ p2 A tdAB p , jyABCD p /q\ 

-n^GB — -tt — 4/1 Kab + -K J^ABCD ■ \p) 

Then the gravity field equation in 6-dimensional is described by the Einstein-Gauss- 
Bonnet equation. 



1 
where 



Gab + oiHab = tti^ab , (9) 



Hab = —7:9abRgb + '^RRab — ^RacR b — 4-R Raced + '^Ra Rbcde ■ (10) 

While the EGB equation Q is rather complicated, it can be shown that in factorizable 
spacetime, the EGB equation can be simplified into a rather illuminating form: after in- 
serting the ansatz (^ into the EGB equations (jHl), the transverse and longitudinal EGB 
equations can be simplified to give 

aRlj,[g]+R[g] = -^J^ , (11) 



{aR[g] + l)RYi] + \R[g] = -^,K + ^1 • (12) 



From those two equations we can see that the main feature of codimension 2 branes in Ein- 
stein gravity is retained in EGB gravity: the scalar curvature of the 4-dimensional spacetime 
is still determined only by the transverse component of the bulk energy-momentum tensor 
from the transverse EGB equation (jllj) : the bulk geometry is then determined by the longitu- 
dinal equation (jl2|) . So the EGB equations are still "factorizable" in factorizable spacetime. 
Thus in the EGB gravity, we still can move the cosmological constant problem completely 
into the bulk. 

Now we have good motivation to proceed to see how the spacetime solutions will be 
modified in EGB gravity. As a first remark, it is interesting to see from Eqs. flllll and 
(|12|) that the Gauss-Bonnet term couples only with the 4-dimensional scalar curvature R[g\. 
Thus, when the 4-dimensional geometry is flat, EGB equations will always reduce to Einstein 
equations (0) and (0). So in this case, the bulk solution is the same as the one given in 
Ref. [6|. Thus what is really interesting is the case when the brane geometry is not fiat. In 
the next section, we will consider de Sitter geometry on the brane. 

III. DE SITTER BRANES IN EINSTEIN-GAUSS-BONNET GRAVITY 

From the recent cosmological observation that our universe is currently accelerating J2l| , 
we are interested in solutions with de Sitter geometry on the brane. So we will consider 
in this section the case that the geometry on the brane is de Sitter, i.e. R[g\^u = ^aQ^iu 
and R[g\ = 4A4, where A4 > is the 4-dimensional cosmological constant. Under those 
assumptions of spacetime geometry, it can be seen from Eqs. (fTT|) and (jT^ that the bulk 
energy-momentum tensor T^^ must be constant along the bulk. 

Let's first discuss the 4-dimensional geometry by Eq. fjllll . Under the assumption of 
maximal symmetric, it can be rewritten as 

^aAl + 4A4 = -^Tf , (13) 

From Eq.()13|). we can find that the 4-dimensional cosmological constant is given in terms 
of the bulk energy-momentum tensor by 

3 



''' 4„ 



-'^'^'-ik^^' 



'6 



(14) 



Thus, the first different feature we encounter in the EGB gravity is that, for any given 
rf^, unless it satisfies T^ = -^^, 



T45, unless it satisfies T^ = ^-^, we will have two solutions of the brane geometry. After 



the brane geometry is determined, the bulk geometry is uniquely determined by the brane 
geometry from Eq. ()12|) . Thus, generally, for any given bulk matter content, there will be two 
different solutions of the EGB equation. This is obviously not a pleasant feature. However, 
we will argue that the "-" branch of the solution is unphysical and should be discarded. 
It can be seen from Eq. (jl4|) that for the "-" branch, the coefficient of -R[7] in Eq. (jl2|) . i.e. 
4aA4 + 1, is always negative; while in Einstein gravity, i.e. a = 0, it is always positive. This 
means that for the "-" branch, the gravity in the transverse dimension is repulsive: positive 
bulk energy density will give rise to negative curvature and only negative tension branes can 
give rise to a positive deficit angle. We think those properties are too exotic so should be 
regarded as unphysical. Thus in the following discussions, we will discard the "-" branch. 
Then let's discuss the bulk geometry from Eq. (jl2j) . which now can be written as 

(4aA4 + l)M^Rb] = -^Tf - 2M^A, + ^S^'\y) . (15) 

|7| 



First, in the case of a vacuum bulk, i.e. T^^ = 0. From the "+" branch of Eq. ()14|l . we 
have A4 = 0, and Eq. p5|) will just reduce to Einstein gravity. Thus the bulk geometry will 
be same as the case discussed in Ref.|l9|. 

Next, let's consider the presence of bulk fields. Following Refs.jy, |7|, |20|, we expect the 



extra dimensions to have the topology of a sphere S^. Thus the 2-dimensional metric '-fab 
will be of the form, 

labdy'^dy^ = aKdO"^ + /^^ sin^ Odip^) , (16) 

where oq is the size of the extra dimensions and /3 is related to the deficit angle 6 by 
5 = 27r(l-/5). 

Transforming the metric (fT^ into the conformal form, 

^abdy'^dy' = ij{r){dr' + r'd^') , (17) 



where ip is given by 



4/5^00 . . 

'^y'^) - r2[(r/ro)^ + (r/ro)-^]2 ' ^ ^ 



and substituting this into Eq.(fT3j). it can be found that ao and P are given by 

, M^n+iaAA 



■\TB - MiA, 



(19) 



^"^ 27rM|(l + 4aA4) ■ ^^°^ 



Equations ()19p , ()2U|) and ()14|) determine the brane and bulk geometry completely. They are 
the main result of this paper. Below we will discuss mainly its application to the scenario 
of Ref.p. Before that, two remarks are in order about those solutions. 

First, a whole new feature of the solution (J2(J|) compared to the Einstein case is that the 
deficit angle in the extra dimensions will now depend on the geometry of the branes. From 
this, we can find an interesting geometric argument in favor of a positive Gauss-Bonnet 
coupling constant. In the case of a negative Gauss-Bonnet coupling constant, the geometry 
will become singular when A4 > — l/(4a). Since we expect A4 to be very large during the 
inflation era, the requirement of a nonsingular geometry forces us to rule out a negative 
Gauss-Bonnet coupling constant. 

Second, the brane geometry in 6-dimensional Einstein-Gauss-Bonnet gravity is also dis- 



cussed in Ref. 



Q 



The authors actually considered only the longitudinal component of the 
EGB equation and concluded that Einstein gravity will restore on the brane. Due to our 
analysis, the brane geometry is determined by the transverse component of the EGB equa- 
tion and while the longitudinal equation looks like Einstein equation, it actually determines 
the bulk geometry after the brane geometry is found by the transverse equation. This can be 



seen more clearly by the expression for the 4-dimensional cosmological constant in Ref. [11 1 
(Eq.(21) in that reference). Actually, Eq.(21) in Ref.|lj| is exactly Eq. ()20j) . from which we 
can see that it actually determines the deficit angle after the 4-dimensional cosmological 
constant is found from Eq. (|Tl|) . 

Now, let's discuss a specific example of the solutions (fnijl . (j^ and ((13)). A lot of the recent 
works on codimension 2 branes are motivated by the exact solution presented by Carroll 
and Guica J6[ which shows explicitly the independence of the 4-dimensional geometry on the 
brane tension. Thus we think it is most important to discuss the corresponding solutions 
in EGB gravity and compare it with that of Ref.p. The solution presented by Carroll and 
Guica assumes a bulk cosmological constant and a magnetic flux, which is described by the 
bulk action p, 

Se = J d^X^l (Im^R - a - IfabF^'A , (21) 

where Mq is the 6-dimensional reduced Planck mass and A is the 6-dimensional vacuum 
energy density. The 2-form field strength takes the form Fab = \^\l\Bo£ab, where Bq is a 
constant and eab is the standard antisymmetric tensor. Other components of Fap vanish 
identically. This model is originally suggested to stabilize the extra dimensions 




The bulk energy-momentum tensor contains contributions from both the bulk cosmolog- 
ical constant and the gauge field, 

T^B = T\b + T^B , (22) 

for which the explicit forms are 



rpX 

-l-AB — ~ 

nB = ABl['"' ' |. (23) 

-lab , 

So we have Tf = -4A - 2Bl and Tf = -2 A + fi^. 

At first, we generalize the fiat brane solution of Ref.jy] to de Sitter brane, which is given 

by 

a- - ^^' (24) 

"^0 " 2A - 3M|A4 ' ^ ' 

MIK, = Ia - \bI . (26) 




It is interesting to note that for the geometry to be nonsingular, from Eq. (j24j) . we must 
have A4 < 2A/(3Mg). However, from Eq.(|^^. this is always satisfied. Thus the de Sitter 
geometry of the brane will never make the bulk geometry singular. 

From Eq. (j26p . we can see that the puzzle of a small 4-dimensional cosmological constant 
is now transformed to the question of explaining a fine-tuning between the 6-dimensional 
vacuum energy and the magnetic fiux, which is a purely bulk problem. Thus in this scenario 
the cosmological constant problem is moved completely into the bulk. Of course, this does 
not solve the cosmological constant problem, but it transforms the nature of the problem in 
an interesting way. At a first glance, it is tempting to appeal to the usual supersymmetry 
argument [^ to set both A and B^ very small, thus avoiding fine-tuning between them. 
However, this cannot work. From Eq.(|24p. we can see that we must require either A or B'^ 
to be of the order Mg so that the size of the extra dimensions can be phenomenologically 
viable. Thus, there is a real fine-tuning problem in the bulk. Currently, we still do not know 
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whether this fine-tuning can be technically natural. Thus, it would be very interesting that 
if in the EGB gravity, we can have a way to release this fine-tuning. We will see below that 
when the Gauss-Bonnet coupling constant is large, this is possible. 

Then, we turn to the discussion of solutions in EGB gravity. From Eqs. (jl9|) . (j^Uj) and 
(J14p . the corresponding solution in EGB gravity is given by 

2 M^il + AaKA 



2A - 3MiA4 - laMiA 



4a2 ' 



(27) 



^"^ 27rM|(l + 4aA4) ' ^^^^ 




3 Ml -') • (2^) 

As a first remark, while we have discussed above, a negative a may result in a singular 
spacetime; for a positive a, while it is not very obvious, it still can be shown that the 
geometry is always nonsingular by an argument that is similar to the Einstein case. 

Then let's discuss the cosmological constant problem in EGB gravity as expressed by 
Eq.()29p. Since we generally have A < Mg and Bq < Afg, so when aM^ < 1, we have 
a(2A - Bl)/Ml < 1. Expanding the RHS of Eq.dlHl) to first order, we can find that 
A4 ~ (2A — Bq)/Mq. Thus, in this case we are actually facing the same fine-tuning as 
in Einstein gravity in order to get a small cosmological constant. This is not a surprise, 
since it is natural for the solution to reduce to the Einstein case when a is small. So 
what is interesting is the case where the Gauss-Bonnet coupling constant is large. Let's 
assume 2A — -Bg ~ Mq, i.e. we do not have a fine-tuning in the bulk, and when aM^ ^ 1, 
i.e. considering the case of a large Gauss-Bonnet coupling constant, from Eq.(|^. we can 
obtain 

A4 - ^ . (30) 



Thus even if we do not have a fine-tuning in the bulk, for a sufficiently large a, we still 
can get a small 4-dimensional cosmological constant. In this case, the current cosmological 
expansion acceleration is actually driven by the 6-dimensional Gauss-Bonnet term, which is 
in some sense similar to the recent model of 1/R gravity proposed by Carroll et al. J2^: the 
current cosmological expansion acceleration is driven by a 1/R term in the 4-dimensional 
gravitational lagrangian. Of course, in order for the A4 to be the order of the observational 



value, a~^ also needs to be fine-tuned to an extremely small value. Thus in the current case, 
we have actually traded the fine-tuning in the bulk to a fine-tuning in the Gauss-Bonnet 
coupling constant. Although this new fine-tuning also seems unnatural now, the cosmological 
constant problem is so hard to solve that it is worth transforming it to a new problem for 
further investigations. Furthermore, this shows the qualitative feature of what will happen 
if we consider higher-derivative gravity in the bulk. Maybe considering more complicated 
higher- derivative gravity theories, such as forth order combinations of the curvature tensor 
can further release the fine-tuning in a more natural way. This deserves further investigating. 
It is worth mentioning that similar fine-tuning problem also happens in the 1/R gravity: 
the coefficient of the 1/R term should also be extremely small to account for the current 
cosmic accelerating expansion [22]. In the 1/R gravity, this is unnatural from an effective 
field point of view and can lead to some inconsistencies when the theory is treated quantum 
mechanically J23|. Now we still do not know whether similar problem will be presenting 
here. 

As a final remark, in Ref.J2J], Navarro considered using a 4-form field in place of the 
2-form field in the action 1)211) . By using Eqs. ()19j) . ()2())1 and (fT^ . it is trivial to generalize 
Navarro's solution to the EGB gravity. 

IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have discussed the gravitational properties of codimension 2 branes in 
Einstein-Gauss-Bonnet gravity and their implications in addressing the cosmological con- 
stant problems. 

Although the current scenario is originally introduced to discuss the cosmological con- 
stant problem, it is also mandatory that cosmological models from String Theory should be 
reconciled with inflation, now a quite well-established ingredient of modern cosmology [22| 
(see, e.g., Ref. 25| for a recent review of braneworld inflation; inflation in 5-dimensional EGB 
gravity is recently discussed in Ref.J26|). When considering inflation model building in the 
present scenario, an observation is that the inflaton must be a bulk fleld. This is in sharp 



contrast to the discussions of the codimension 
assumes the inflaton to be conflned on the brane 



(lase, where most of the inflation model 
. The reason for this is simple. Current 



observation of the CMB power spectrum tells us that during inflation, the energy density 
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of inflaton should be almost constant [27]. Thus, if the inflaton is a field confined on the 
brane, then during inflation it will behave just like the brane tension. So the above analysis 
tells us that it cannot affect the 4- dimensional geometry. On the other hand, if the inflaton 
is a bulk field, then its effects during inflation are just equivalent to a renormalization of the 
6-dimensional cosmological constant A. 

So in the Einstein gravity case, the Hubble parameter during inflation H^ = A4/3 will 
be given from Eq. lj^UI) by 

H' ~ V/M^ , (31) 

where V is the potential of the bulk inflaton fleld. Thus the energy scale of the potential 
would be of order {H/MqY'^Mq during inflation. In the original model of large extra dimen- 
sions yl, in order to address the gauge hierarchy problem, the 6-dimensional reduced Planck 
mass is assumed at most a few orders higher than the supersymmetry breaking scale which 
is of order 1 TeV. On the other hand, current CMB data prefers a high inflation scale which 
is at most several orders of magnitude smaller than the GUT scale ~ 10^^ GeV [27] . Thus, 
the potential V during inflation is necessarily larger than Mg. If the inflaton is a brane fleld, 
there is nothing unnatural here. But as we have commented above, inflaton must be a bulk 
fleld now. So it is very unnatural for a bulk fleld to have an energy scale larger than the bulk 
Planck mass. Therefore, implementing successful inflation scenario encounters fundamental 
difficulties in codimension 2 brane scenarios |3ll |. 

The situation is worse in EGB gravity. From Eq. ()29j) . the Hubble parameter during 
infiation will be given by Eq. ()3Hl when M^a <^ 1 and it reduces to the Einstein case. When 
M^a ^ 1, from Eq. ()30|) . the Hubble parameter will be given by 

H^ ^ Ml^ (32) 



Thus the energy scale of the potential would be of order (if/M6)^/^(aM|)^/^M6 during 
infiation. Then a higher potential energy is needed compared with the Einstein case ()31|) to 
implement the infiation. This makes the problem we discussed above more severe. 

Faced with the above problem, it is worth considering other mechanism of driving an 
infiation on the brane rather than a bulk scalar field. A seemingly promising candidate is 
the R^ infiationary model of Starobinsky [2^. However, in order to avoid the above problem, 
we assume that the R^ term is only induced on the brane, like the case of induced gravity 
model given by Dvali et al. J29|. More concretely, we may consider adding to the bulk 
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lagrangian (PT|) an induced R^ term, 

Sinduced = / d^x^/\g\aR[g]^ , (33) 



where a will be of order M^'^ J28| . It is worth commenting that such a term may be induced 
by quantum effects of conformal fields on the brane, and R^ inflation on codimension 1 
braneworld has been discussed in Ref.J30|. This and other possibilities to handle the inflation 
model building problems deserve further investigation. 
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